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There are well-known protocols for performing CNOT quantum logic with qubits coupled by 
particular high-symmetry (Ising or Heisenberg) interactions. However, many architectures being 
considered for quantum computation involve qubits or qubits and resonators coupled by more com- 
plicated and less symmetric interactions. Here we consider a widely applicable model of weakly but 
otherwise arbitrarily coupled two-level systems, and use quantum gate design techniques to derive a 
simple and intuitive CNOT construction. Useful variations and extensions of the solution are given 
for common special cases. 
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Experimental realizations of gate-based quantum com- 
putation require the accurate implementation of uni- 
versal two-qubit operations such as the controUed-NOT 
(CNOT) quantum logic gate [H. Finding the best way 
to achieve this for a specific experimental architecture is 
a principal goal of what we refer to as quantum gate de- 
sign. The CNOT problem can be informally stated as 
follows: Specify the dimension N of the relevant Hilbert 
space, and a Hamiltonian 



(1) 



with some experimental control over K parameters 
■ ■ . ,S,K- How should the control parameters be var- 
ied to generate CNOT logic in the computational ba- 
sis {|00), |01), |10), |11)}? For a closed system this is a 
control problem in the unitary group U(A^). N is not 
necessarily equal to 4 because the Hamiltonian might 
include auxiliary non-qubit states (not in the computa- 
tional basis) that help implement the logic. For exam- 
ple, an effective strategy (see, for example, Strauch et 
al. [1|) is to use an anticrossing of the |11) state with 
a non-computational state |aux) to generate a 27r rota- 
tion in the two-dimensional subspace {|11), |aux)}. This 
implements the gate CZ = diag(l, 1,1,-1) in the com- 
putational basis, out of which a CNOT can be made 
by pre- and post-application of single-qubit Hadamards. 
Another important example is Cirac and Zoller's use of 
vibrational modes to mediate quantum logic between the 
internal qubit states of trapped ions 13]. 

In the familiar U(4) case, the Hamiltonian ([1]) can be 
written in terms of Pauli matrices and their tensor prod- 
ucts. The resulting coupled- qubit model usually allows 
control of some of the single-qubit operators — enough to 
perform arbitrary SU(2) rotations on each qubit — and 
possibly of the qubit-qubit coupling. For certain com- 
monly occurring forms of the qubit-qubit interaction, in- 
cluding the highly symmetric cases of Ising-like crla2 in- 
teraction j| and Heisenberg-like cri ■ ct2 interaction [3], 
effective protocols for implementing CNOT gates have 
been established. However, many architectures being 



considered for quantum computation involve qubits or 
qubits and resonators coupled by more complicated and 
less symmetric interactions, or would be more accurately 
modeled as such. 

Here we investigate the general problem of weakly 
but otherwise arbitrarily coupled qubits, and use per- 
turbation theory combined with other quantum gate de- 
sign techniques to derive a simple and widely applicable 
CNOT pulse construction. Useful variations and exten- 
sions of the basic solution are given for common special 
cases, and the intuitive geometric picture we employ (re- 
lated to the Weyl chamber description used by Zhang et 
al. [5]) will be useful elsewhere in the design of quan- 
tum logic gates. We assume unitary evolution, which is 
sensible given the generality of our result and the wide 
variation in experimental coherence times. 

Zhang and Whaley Q have addressed the problem of 
two-qubit gate construction using similar methods ap- 
plied to a variety of coupled-qubit models, but focused 
on steering with continuous rf control as opposed to the 
short pulses considered below. One of us has recently 
investigated the implementation of CNOT gates using 
constant rf driving [3, Q and moderately detuned qubits 
[H , providing constructions complementary to those pre- 
sented here. Time-optimal and other direct quantum 
control approaches are especially useful for strongly cou- 
pled and/or strongly driven qubits, or to optimize per- 
formance in the presence of specific decohering and/or 
noisy environments [III . [13 . [ist , but early quantum 
logic demonstrations might best be accomplished using 
the simple perturbative protocol described here. 

Weakly coupled qubits. In a wide variety of physical 
systems being considered for quantum computation, the 
Hamiltonian for a pair of coupled qubits can be written 
(suppressing K) as 

^ = E (-f'Tf + ^^.cos(e^^ + 0,)<) 
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with J^jy a 3x3 real-valued tensor (possibly adjustable). 
The Hamiltonian ^ is written in the basis of eigenstates 
(|0) and |1)) of uncoupled qubits with energy level spac- 
ings Ei, and the parameters and f2i (with <^ e^) are 
assumed to be experimentally controllable. More gen- 
eral single-qubit control is often available but will not be 
needed here. Our principal assumption is that of weak 
coupling: The magnitude of the J^^, are assumed to be 
small compared with the Ci. 

Two-qubit logic gates will be implemented by combin- 
ing certain entangling operations, performed with tuned 
(ei = £2) qubits, together with single-qubit operations 
performed with detuned or decoupled qubits [14|. In a 
frame rotating with the tune d q ubits, the Hamiltonian 
([2]) reduces approximately to 
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where 



n 



Here 



Jxx ~^ Jyy Q^Yi'di j' '^^y '^yx 



(5) 



In the computational basis. 
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where 7 = 2{J + iJ') = \j\ e**^. To obtain ^ we have as- 
sumed that the J^i^ and Q are small compared with the 
qubit frequency and have neglected the resulting rapidly 
oscillating terms with vanishing time-averages (the usual 
rotating- wave approximation). Although 9 coupling con- 
stants are present in ([2]), only 3 parameters appear in 
Ti., making a general analysis possible. The terms in 
(III) multiplying J and Jzz are symmetric under qubit- 
label exchange, whereas the J' term is antisymmetric and 
therefore vanishes when the physical qubits in question 
(and their operating biases) are identical. Furthermore, 
in the common case of J' = (which must occur when 
the qubits are identical but can also occur when they are 
not), H commutes with itself at different times when J 
and Jzz are time dependent, leading to additional flexi- 
bility (in the form of "area" theorems) for pulse design 
that we will use below. We emphasize that 7i is a univer- 
sal Hamiltonian, applying to any pair of tuned, weakly 
coupled qubits. Coupled-qubit models with nondiagonal 
single-qubit drift terms can be put in the form ^ after 
transformation to the uncoupled eigenstate basis. 
Cartan decomposition. The trajectory in U(4) that 



traces out during Schrodinger evolution (T is the stan- 
dard time-ordering operator) can be viewed by factoring 
out local (single-qubit) rotations u G SU(2) (giSU(2). A 
convenient way to achieve this is to use the fact that any 
element of U(4) can be written as 



with 



A{x,y,z) 
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for some local rotations Uprc and Wpost, real- valued coor- 
dinates (angles) x, y, and z, and global phase (j). This 
formula can be derived by using a Cartan decomposition 
of the Lie algebra su{4) [I^, IQ, 17 1. The central com- 
ponent A has the geometrical structure of a 3-torus with 
period 2tt and characterizes the nonlocal or entangling 
part of U. By performing the decomposition ([8]) at each 
time t and forming the vector r = (x, y, z), we can view 
the evolution of the nonlocal part of C/ as a trajectory r(t) 
through the three-dimensional space of entanglers ([9]). A 
special property of ([9|) is that the generators cr^ , 
(jy® , and cr^® all commute (they form the abelian 
Cartan subalgebra). The minus sign introduced into the 
exponent of ^ simplifies the analysis in the common 
special case of J' = 0. 

The decomposition ([8]) into an entangler A, local rota- 
tions Upro and Upostj and phase factor e"^ is not unique. 
This means that the trajectory f{t) corresponding to 
some actual physical evolution is not unique. But the 
different options for A at each time t are evidently lo- 
cally equivalent (differing by pre- and post-application of 
local rotations and a multiplicative phase factor). Fur- 
thermore, in the common special case of J' = 0, a particu- 
larly natural continuous solution [given in (jlip below] can 
always be chosen which has the simplifying property that 
the local rotations and phase factor are equal to the iden- 
tity along the entire trajectory: The local rotation and 
global phase angles vanish. We note that the usefulness 
of the decomposition ((8]) goes far beyond its somewhat 
technical role here: (i) In architectures where local opera- 
tions can be performed quickly and accurately (they are 
"free"), the decomposition allows one to focus directly 
on the remaining nonlocal part; (ii) The local rotations 
associated with successive gates can often be combined; 
And (iii), some of the experimental error incurred when 
implementing an entangler — the component that doesn't 
change the equivalence class — can be corrected by mod- 
ifying the u's. 

The concepts of local equivalence and local equivalence 
classes have wide application in gate design. Makhlin 
has constructed an explicit formula for 3 quantities that 
can be used to test for local equivalence. The CNOT 
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FIG. 1: (color online) Three-dimensional space of entanglers 
A, showing the six members ^(±^,0,0), ^(0,±-|,0), and 
j4(0, 0, ±-|) of the CNOT equivalence class (solid green cir- 
cles) closest to the identity (open red circle). The Schrodinger 
evolution resulting from Q with fixed positive J and Jzz is 
indicated by the black trajectory, interrupted by the applica- 
tion of a fast n pulse. The purple square at r = (f , fiO) is 
locally equivalent to the CNOT x SWAP and SWAP x CNOT 
gates, and the blue star at r = (f , f, f ) is locally equivalent 
to the SWAP. 



lustrated in Fig. [T] The fact that the coordinates 
([TT|) of the generated entangler depend only on time- 
integrals of the coupling constants indicates a type of 
robustness and flexibihty of the associated experimen- 
tal pulse sequence, analogous to the area theorem for 
single-qubit rotations within the rotating-wave approx- 
imation. The closest CNOT-class entanglers are at 
f = (±f,0,0), (0,±f,0), and (0,0,±f), which can- 
not be reached with one application of Ti unless J van- 
ishes. In this Ising case, A{0,0,±j) is obtained after 
/ Jzz dt = ±^ (mod27r). Generating one of these entan- 
glers corresponds to generating a particular member of 
the CNOT equivalence class; either one is sufficient. A 
possible Ising pulse sequence (executed from right to left) 
is CNOT - eT^f H2i?.(Tf)ii?z(Ti)2^(0,0^f)H2, 
where H = i Rx{t^) Ry{'§) is a Hadamard gate [2l|. 

Another important special case occurs when J^z van- 
ishes, often called an XY interaction Here one 
can follow the general "two-shot" protocol detailed be- 
low to generate the canonical CNOT gate (|10[) or, alter- 
natively, one can generate A{±j,±j,0) in a single shot, 
which is locally equivalent to both CNOT x SWAP and 
SWAP X CNOT, where 
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has Makhlin invariants Gi — and G2 = 1 (Gi is gener- 
ally complex). Two members U and U' of U(4) are locally 
equivalent if and only if their Makhlin invariants are iden- 
tical, in which case we write U ^ U' . When restricted to 
a certain (nearly) tetrahedral region — a Weyl chamber — 
the angles x, y, and z are in one-to-one correspondence 
with the Makhlin invariants, leading to a unique r and 
a beautiful geometric description of the local equivalence 
classes of U(4) [5]. For our purposes, however, it will be 
convenient to work in the full toroidal space of entanglers 
and not restrict r to a Weyl chamber [20| . 

CNOT construction when J' — Q. First we consider 
the common special case of Hamiltonian ^ with J' = 0, 
which includes the case of identical qubits. Assuming 
tuned qubits and no rf drive, the evolution ([7|) simplifies 
ioU = A, with 



Jdt, Jdt, Jzzdt , (11) 



which follows a curve in the vertical plane x = y. 
The trajectory for the case of fixed J and Jzz is il- 



The gates CNOT x SWAP and SWAP x CNOT are as 
effective as (jlOp in the sense that any quantum circuit 
written in terms of CNOTs can be immediately rewritten 
in terms of the swapped versions with no overhead [23j . 
An example pulse sequence is 



SWAP X CNOT === ±z ± f ) ^ i?4f ) ^] [i?, (f ) ^ 



X Ry 



J^4-f)2^(±f'±f'0) 



2)2- 



(13) 



Operations grouped together in square brackets can be 
performed simultaneously. 

Although CNOT-class entanglers farther from the ori- 
gin (and not shown in Fig. [T]) can be reached in one shot 
for special values of Jzzj J a faster and generally appli- 
cable protocol is to interrupt the evolution with a fast 
refocusing tt pulse applied to either qubit. A pair of such 
pulses enclosing an interval of tuned qubit evolution 



i?,(-^)e-»/«''*i?,(7r) 



(14) 



can be viewed as transforming the interaction Hamilto- 
nian during that interval to (note sign changes) 



1 { 



(15) 



causing the reflection illustrated in Fig. [T] and allowing 
the evolution to reach any entangler on the positive x 
axis (or negative axis for J < 0) . Ry{T^) would cause a 
reflection toward the y axis. 
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Entanglers on the x, y and z axes are locally equivalent 
to each other and to the controUed-phase gate 
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^A(|,0,0)~A(0,|, 0)^^(0,0,. 



(16) 

The entanglers at r = (if, 0, 0), shown in Fig.[l]as solid 
green circles, are thus locally equivalent to the CZ gate, 
and hence to the CNOT. The identity 

i?4-7r)^A(±f,±|,z) i?47r)^=A(±f,Tf,-^), 

with z arbitrary, allows us to reach 

A(±f,Tf,-^)A(±f,±f,^) =^(±f,0,0) (17) 

after two entangling intervals, out of which a CNOT can 
be constructed according to 

CNOT = eT'ii?,(-f)^[i?4Tf)ii?.(Tf)2] 

X A(±f,0,0)i?,(f)^. (18) 

Arbitrary J'. Here we assume Hamiltonian Q with 
fixed, time-independent values of J, J^^, and J' (ex- 
cluding the pure erf erf Ising case). When J' / there 
are terms in the Hamiltonian that are not in the Cartan 
subalgebra and that break the symmetry under qubit ex- 
change. Such terms can be eliminated by performing a z 
rotation on the second qubit by an angle ip = arg( J-l-i J'), 



+ Jzzcjlal, (19) 
allowing us to reach the CNOT-class entangler 
A(f,0,0) = R.{-^)^R,{~t:)^ e-'«^* 



-iHAt 



(20) 



Here e ^* represents the action of bringing the qubits 
into resonance for a time 



At; 



(21) 



The complete pulse sequence in this case can be written 
as 



Brown, Sue Coppersmith, Andrew Sornborger, Matthias 
Steffen, and Robert Varley for useful discussions. 



[2] 

[3] 
[4] 

[5] 

[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 



M. A. Neilsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, England, 2000). 

F. W. Strauch, P. R. Johnson, A. J. Dragt, C. J. Lobb, 
J. R. Anderson, and F. C. Wellstood, Phys. Rev. Lett. 
91, 167005 (2003). 

J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 
D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

J. Zhang, J. Vala, S. Sastry, and K. B. Whaley, Phys. 
Rev. A 67, 042313 (2003). 

J. Zhang and K. B. Whaley, Phys. Rev. A 71, 052317 
(2005). 

A. Gahautdinov, Phys. Rev. A 75, 052303 (2007). 

A. Gahautdinov, J. Math. Phys. 48, 112105 (2007). 

A. Gahautdinov, Phys. Rev. A 79, 042316 (2009). 

N. Khaneja, R. Brockett, and S. J. Glaser, Phys. Rev. A 

63, 032308 (2001). 

P. Zanardi and S. Lloyd, Phys. Rev. A 69, 022313 (2004). 
I. A. Grigorenko and D. V. Khveshchenko, Phys. Rev. 
Lett. 95, 110501 (2005). 

A. Sporl, T. Schulte-Herbriiggen, S. J. Glaser, 
V. Bergholm, M. J. Storcz, J. Ferber, and F. K. Wil- 
helm, Phys. Rev. A 75, 012302 (2007). 
When Jzz ^ 0, the interaction between detuned qubits 
is not suppressed, and it is preferable to switch off the 
coupling during single-qubit operations in this case. 
Transformation to a frame rotating with qubit frequen- 
cies ei and e2 is defined according to \'ip) xl^), 
where x = Yli {J '^^) ■ Under this transformation, 
< ^ X = dt) < , where 



7^.(a) 



cos a sin a i 
— sin a cos a 
1; 



[16] 
[17] 
[18] 
[19] 



[20] 

[21] 
[22] 



CNOT = _ 1)^ _ 1)^] R^{^v)^ 

X 1 e'^'' (^) 1 e-*'' ^* i?. (^) 2 i?, (f ) , . (22) [23] 

In conclusion, we have shown how to implement the 
CNOT quantum logic gate with weakly but otherwise 
arbitrarily coupled qubits. This work was supported by 
lARPA under grant W911NF-08- 1-0336 and by the NSF 
under grant CMS-0404031. It is a pleasure to thank Ken 



is an SO(3) rotation matrix. 

N. Khaneja and S. J. Glaser, Chem. Phys. 267, 11 (2001). 
B. Kraus and J. I. Cirac, Phys. Rev. A 63, 062309 (2001). 
Y. Makhlin, Quant. Info. Proc. 1, 243 (2002). 
The CNOT operator defined in (|10|) specifies qubit 1 to 
be the control qubit. The locally equivalent alternative 
controlling the second qubit is SWAP x CNOT x SWAP. 
We also adopt a slightly different notation than that of 
Ref. 5], writing c as — 2r. 
R^{0)i = e~''*''^'°"i is a, rotation on qubit i. 
Pulse sequences for the XY case were previously derived 
in unpublished notes by Matthias Steffen. 
These gates are also equivalent double CNOTs, consist- 
ing of a pair targeting each qubit in succession (see [Q): 

SWAP X CNOT = CNOT X (SWAP CNOT SWAP) 
CNOT X SWAP = (SWAP CNOT SWAP) X CNOT. 



